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' Abstract. We perform an exact and asymptotic analysis of the model of n vicious 

, walkers interacting with a wall via contact potentials, a model introduced by Brak, 

' Essam and Owczarek. More specifically, we study the partition function of watermelon 

, configurations which start on the wall, but may end at arbitrary height, and their mean 

number of contacts with the wall. We improve and extend the earlier (partially non- 
rigorous) results by Brak, Essam and Owczarek, providing new exact results, and more 
precise and more general asymptotic results, in particular full asymptotic expansions 
. , for the partition function and the mean number of contacts. Furthermore, we relate 

' this circle of problems to earlier results in the combinatorial and statistical literature. 

(N 
m 
\o 

x/^ '. 1. Introduction. The problem of vicious walkers was introduced by Fisher [17], 

O I who also gave a number of physical applications of the model, such as, for example, 
^ ' to modelling wetting and melting. The general model is one of n random walkers on 
' a (i-dimensional lattice who at regular time intervals simultaneously take one step in 
the direction of one of the allowed lattice vectors such that at no time two walkers 
occupy the same lattice site. 
^ ; Numerous papers have been written on the subject since then. Most of them 

■ analyse the model of vicious walkers in a continuum limit (such as for example [17, 



20, 21, 22]). It has been realized only recently that in fact there are many interesting 
cases in which even exact results in form of nice closed product formulas are available 
(see for example [1, 7, 8, 9, 15, 29, 37, 41, 45]), and that asymptotic analysis can be 
performed directly on the model, without taking recourse to continuum limits, thus 
obtaining more precise estimates (see for example [8, 15, 37, 41, 45]). 
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In this paper we consider n vicious walkers (the "vicious" constraint demanding 
that at no time two walkers occupy the same site) in the plane integer lattice with 
allowed steps of the form (1,1) (up-step) and (1,-1) (down-step), which in addition 
do not run below the x-axis (the wait). Given integer vectors a = (ai, 02, . . . , a^) 
and e = (ei, 62, . . . , e^), where all a^'s are of the same parity and all e^'s are of the 
same parity, we study the partition function z[^\ai — > e; k) := where the 

sum is over all families P = (Pi, P25 ■ ■ • , Pn) of n vicious walkers as above, the i- 
th walker Pi starting at (0, a^) and ending at (t, e^), i = 1,2, ... ,n, and with c(P) 
denoting the total number of contacts of the walkers with the wall. (Because of the 
"vicious" constraint, it is only the lowest of the walkers who can have contacts with 
the wall.) An example of such a family of vicious walkers with n = 4, a = (0, 2, 4, 6), 
e = (0,2,4,6), t = 8 is shown in Figure l.a, whereas Figure l.b shows a family of 
vicious walkers with n = 5, a = (0, 2, 4, 6, 8), e = (2, 4, 6, 8, 10), t = 12. The number 
of contacts with the wall (the x-axis) is 3 in Figure l.a, and it is 4 in Figure l.b. 




The special topology of vicious walkers that we shall be mostly concerned with in 
this paper is watermelon configurations {watermelons for short), which are families 
of vicious walkers in which the ^-coordinates of neighbouring starting points differ 
by 2, the same being true for the end points. The y-coordinate of the end point of 
the lowest walker is called the deviation of the watermelon configuration. Thus, both 
vicious walker families in Figure 1 are in fact watermelon configurations, the one in 
Figure l.a has deviation 0, the one in Figure l.b has deviation 2. 
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The above described vicious walker model with contact interaction has been in- 
troduced by Owczarek, Essam and Brak in [41]. They obtained a dctcrminantal 
formula for the partition function, and undertook a (partially non-rigorous^) study 
of its asymptotic properties as the length t of the walks tends to infinity in the case 
of watermelon configurations starting and ending on the wall. (By definition, these 
arise if a = e = (0, 2, 4, . . . , 2n — 2).) This study was crucially based on recurrence 
relations for the partition function. It revealed a critical value of the parameter k at 
K = 2, and they also provided a scaling analysis around this critical value. However, 
the use of recurrence relations intrinsically makes it only possible to obtain the order 
of magnitude of the partition function, but not the multiplicative constant, not to 
mention any error terms. This paper was followed by the articles [6, 7] by Brak and 
Essam, in which exact results for the partition function for watermelons which start 
on the wall and end at some fixed deviation y not necessarily equal to are obtained. 
The latter three papers do on the other hand not contain any asymptotic results for 
this more general situation. 

The purpose of our paper is three-fold. First, we relate this circle of problems 
to earlier results in the combinatorial and statistical literature. In particular, we 
show that the determinantal formulas in [9] and [41] follow directly from a (now) 
classical result on non-intersecting path families in acyclic graphs due to Lindstrom 
[38], which was rediscovered by Gessel and Viennot [24, 25]. (In fact, all determinantal 
formulas for vicious walkers follow from that result.) Furthermore, we outline that 
the results in [8] and [41] on the partition function of a single walker are well-known 
in nonparametric statistics. Second, we derive a new exact formula for the partition 
function 

Zi") {y; k) := Z^^^ ((0, 2, . . . , 2n - 2) ^ (y, y + 2, . . . , y + 2n - 2); «) (1.1) 

for watermelon configurations which start on the wall and end at deviation y (see 

Theorem 8), which expresses zj:^\y; k) in terms of a double sum. This is the central 
result of our paper, from which all other results are derived. Not only does it allow 
us to rederive, in a uniform manner, all previously obtained exact results [6, 7, 8, 

11, 37, 41] on the partition function Z^^\y; k) in the literature (see Corollary 5, 
Theorem 9, Corollary 10, Theorem 11), it also enables us to extend these in two cases 
to significantly larger domains of the parameter k (sec Corollary 5 and Theorem 11). 
Third, we improve and extend the asymptotic results by Owczarek, Essam and Brak 
[41] . Not only do we show how to rigorously find the asymptotic form of the partition 
function for watermelon configurations which start on the wall and end at an arbitrary 
(fixed, or non-fixed) deviation as the length of the walks tends to infinity, our approach 
allows to even derive full asymptotic expansions. In particular, our results confirm 
the phase transition ai k, — 2 predicted by Owczarek, Essam and Brak. Moreover, 
we provide new exact and asymptotic results for the (normalized) mean number of 
contacts of the watermelons with the wall. These results show that, for the length of 



^ Their approach could probably be made rigorous by applying the Birkhoff-Trjitzinsky theory 
[2, 3] of determining the asymptotics of solutions to difference equations, surveyed in [53]. 
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the walks being large, the normalized mean number of contacts is (asymptotically) 
proportional to a constant if k is less than the critical value 2, it is proportional to 
the square root of the length of the walks if k is equal to the critical value, and it is 
proportional to the length of the walks if k is greater than 2. Again, all these results 
are made by possible by our double sum formula for the partition function given in 
Theorem 8. 

The techniques that we employ to prove our results are (1) combinatorial: path 
manipulation, Lindstrom-Gessel-Viennot theorem on non-intersecting lattice paths, 
tableau combinatorics (in particular: jeu de taquin); (2) algebraic- manipulatory: hy- 
pergeometric series identities;^ and for the asymptotic calculations (3) analytic: sin- 
gularity analysis. In particular, our approach clearly shows how the various results 
that have been obtained earlier are connected. 

Our paper is organised as follows. In the next section we address the analysis of a 
single walker with wall interaction, and relate it to the statistical literature, in par- 
ticular to papers by Engelberg [14] and Mohanty [39]. Then, in Section 3, we address 
the case of several walkers. We recall the Lindstrom-Gessel-Viennot theorem and 
demonstrate how it directly implies the determinantal formula [9] by Brak, Essam 
and Owczarek. In addition, in Proposition 3, we also deduce a slightly different de- 
terminantal formula for the partition function for watermelon configurations which 
start on the wall, which will be more convenient for our subsequent computations. In 
Section 4 we restrict our attention to the exact enumeration of watermelon configu- 
rations which start and end on the wall. We provide a new, short proof for a result 
by Brak and Essam [6, Theorem 6], which expresses the partition function in form of 
a single sum. This proof shows in particular that guessing the result is sufficient to 
prove it, the details being filled in by the computer (see Theorem 4 and its proof). 
Moreover, we reprove and extend an alternative expression for the partition function 
found earlier by Owczarek, Essam and Brak [41, Eq. (4.65)] (see Corollary 5). Sec- 
tion 5 begins the analysis of watermelon configurations which start on the wall and 
end at some fixed deviation y not necessarily equal to 0. The main purpose of the 
section is to find a manageable expression for the number of these watermelons with 
a fixed number of contacts with the wall. This is accomplished by showing that these 
are equinumerous with another set of vicious walker families (see Proposition 6), a 
result that has been previously obtained by Brak and Essam [6, Cor. 1]. Our proof 
however is bijective, being based on a modified jeu de taquin, thus solving a problem 
posed in [6] . (We remark that since the first version of the present article was written, 
Rubey [46] has found a completely different bijection, which, in fact, works in a much 
more general setting.) This result is then used in Section 6 to find an exact expression 
in form of a double sum for the partition function zl:^\y; k) for watermelons which 
start on the wall and end at some arbitrary fixed deviation (see Theorem 8), which 
is new and the central result of this paper. By the use of hypergeometric summation 
and transformation formulas we are able to rederive, respectively extend, previously 



■^All the hypergeometric calculations in this paper were carried out using the author's Mathe- 
matica package HYP, which is designed for a convenient handling of hypergeometric series, and is 



available from tittp : //www .mat .univie . ac . at/~kratt 
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obtained alternative expressions due to Brak and Essam (private communication) (see 
Theorems 9 and 11). By the same techniques we are also able to rederive earlier closed 
form results by the author, Guttmann and Viennot [37] in the case that k = 1, and 
by Owczarek, Essam and Brak [41] and Ciucu and the author [11] in the case that 
K = 2, see Corollary 10. Section 7 is devoted to the asymptotic analysis of the parti- 
tion function Z^^\y; k). Our starting point is the double sum formula for z\^\y^ k) 
provided by Theorem 8, which, when combined with the technique of singularity anal- 
ysis, yields full asymptotic expansions for zj:"'\y; k) as the length t of the walks tends 
to infinity (sec Theorem 12). Thus we confirm the predictions by Owczarek, Essam 
and Brak [41], making them more precise at the same time. The final Section 8 is 
devoted to the study of the mean number of contacts. In the same special cases as 
above, namely for k = 1 and k = 2, we are able to obtain simple closed formulas, see 
Theorems 13 and 14. Singularity analysis again allows us to obtain full asymptotic 
expansions for the mean number of contacts, which is the contents of Theorem 15. 

2. The combinatorics of a single walker. In this section wc consider a single 
walker with wall interaction. What we aim at is the computation of the partition 
function k^^^\ where the sum is over all walks P from (0, a) to (t, e) which never 
run below the x-axis (the wall), and where c(P) denotes the number of contacts 
of the walk with the wall. Slightly deviating from the notation introduced in the 
Introduction, we denote this partition function by Z^^\a — > e; k). 

As we are going to develop, this type of question is classical in the literature on 
nonparametric statistics. We summarize what is known in the proposition below. It 
gives an explicit formula for the partition function zj:^\a — > e;^), and, as a simple 
corollary, an algebraic expression for the corresponding generating function. This is 
an old result due to Engelberg [14, Cor. 3.2] (for the case where the starting point is 
on the a;-axis) and Mohanty [39, Cor. 1 (iv)] (for the general case). It has recently 
been rediscovered by Brak, Essam and Owczarek in [8, Eq. (3.11)]. There is an elegant 
combinatorial proof, the origin of which is difficult to track down (see for example 
[31, proof of Lemma 1] for an occurrence of that idea^). Since it is apparently less 
well-known than it should be, we reproduce it below. 

Proposition 1. Let a and e be non-negative integers with a + e = t mod 2. Then 
the partition function zj:^\a — > e; n) of a single vicious walker starting at (0, a) and 
ending at {t, e) and never running below the x-axis is given by 



^The original inspiration for that idea may be due to Csaki and Vincze [13, p. 100]. There it 
appears in a shghtly modified form, and is used there to solve a slightly modified problem, namely 
to find the number of walks which cross the a;-axis a given number of times. 



{t + a- e)/2 




{t-a- e)/2 
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The generating function, the coefficients of which are these partition functions, is 
equal to 



J2 z['\a^e;K)z' = 

t=a+e (mod 2) 



1 /I - Vl - 



2z 



1 / 1 - VI - 4^2 " 



e+a 



VI - 4^2 



2z 



1- 



2k 
2-k; 



VI - 4^2 

1 + 2^vr^^ 



(2.2) 



Proof. Let for the moment a and e be at least 1. The coefficient of k° in '{a 
e; k) counts the walks from (0, a) to (t, e) which do not touch the x-axis. As is classical, 
this number can be obtained from the reflection principle (see e.g. [12, p. 22]). It is 
equal to 

{{t + a- e)/2) ~{{t-a- e)/2) ' ^^'^^ 

which explains the first line in (2.1). 

On the other hand, for £ > 1 the coefficient of z^ counts the walks from (0, a) to 
(t, e) which do not cross the a;-axis and touch it exactly £ times. Let us consider such 
a walk, see Figure 2 for an example in which t = 17, a = 3, e = 2,£ = 3. 




Such a walk is now transformed into a walk from (0, a) to {t — £,e + i) by deleting 
all the steps immediately preceding a touching point from the original walk and gluing 
the walk pieces together. In our example in Figure 2, these steps are indicated by 
thick line segments. Figure 3 shows the result after deletion of these steps. The circles 
should be ignored at this point. 

We claim that this mapping is indeed a bijection between walks from (0, a) to (t, e) 
which do not cross the x-axis and touch it exactly i times and walks from (0, a) to 
{t — £, e + £) which do not run below the horizontal line y — 1 but touch y — 1 at 
least once. For establishing the claim we have to explain how we can reconstruct the 
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After deletion of steps 
Figure 3 



original walk from one of the walks of the second type. In order to accomplish this, 
we just have to find the points in the walk where steps had been deleted. Indeed, 
as is straight-forward to see, for each J, 1 < J < ^, the right- most point on the walk 
which is on the horizontal line y — j is such a point, and these are in fact all of them. 
In our example in Figure 3 these points are circled. 

Hence, the number that we want to find is equal to the number of walks from (0, a) 
to {t — e + £) which do not run below the line y — 1 minus the number of walks from 
(0, a) to (t — £,e + I) which do not run below the line y = 2. By another application 
of the reflection principle, this is 

t-£ \ _ / t-£ 

{t + a-e-2£)/2) ~ \{t - a - e - 2£)/2 

t-£ \ _f t-£ 
{t + a-e-2£)/2) ~ \{t - a - e - 2£ + 2)/2 

which, after cancellation, gives exactly the expression which appears as the coefficient 
of k\ £ > 1, in (2.1). 

The above considerations establish (2.1) for a, e > 1. If a = and e > 1, then 
we argue that the partition function zl:^\o e; k) is equal to KZj:]}^{l — > e; k) 
because the first step in a walk that starts at (0, 0) must necessarily be an up-step, the 
additional factor k, taking into account the touching of the original walk at (0, 0). Now 
we may apply formula (2.1) with a = 1. It is just a matter of simple manipulations to 
convert the obtained expression to (2.1) with a = 0. If a > 1 and e = 0, respectively 
if a = e = 0, one argues similarly. We leave the details to the reader. 

In order to establish the second assertion of the proposition, we must compute the 
generating function 



E 



t=a+e (mod 2) 



(t + 0-6)72) ((t-a-e)/2j ' ^^""^^ 



on the one hand, and 

t=a+e (mod 2) 

on the other. 

We may concentrate on the case where e > a, because it is combinatoriaUy obvious 
that an interchange of a and e gives the same result. Therefore, from now on we 
assume e > a. 

In order to compute the series (2.4), we write it in a telescoping form as 



§ § {{{t + a-e-2j)/2) {{t + a - e - 2j - 2)/2 

t=a+e (mod 2) 



j=0 r>0 



where in the next-to-last line we performed the index transformation t — > 2r+2j+e—a. 

Now we appeal to the well-known fact (see e.g. [40, displayed line before (1.21) 
with = 1], in combination with [40, (1.19)]) that 

E^("r)^-E((":") - (!:r))^-^«"^^ 

where C{x) is the generating function for the Catalan numbers, 

»«=gri,e)-£(P;)-(.i'.))-=^- <■■« 

If we use this in (2.6), then we find that the series (2.4) is equal to 

E 2j+e-a^f 2\2j+e-a+l _ ^ ) ^ ^ \^ ) 

^ ' - l-z2C(z2)2 

/ / , \ e — a / , \ e+a^ 

/ / 1 - Vl - 4^2 \ / 1 - Vl - 4^2 



1 

Vl - 4^2 



V 



2z ] \ 2z 



(2.9) 



the last line being due to the easily verified fact that 

l-xC{xf = C(a;)Vl -4a; (2.10) 
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For the computation of the series (2.5) we proceed in a similar fashion. Here, 
we first interchange summations and then replace the index t of the (now) inner 
summation by 2r + 2£ + a + e — 2 to transform the series to 



£>1 r>0 ^ ^ 



2r + £ + a + e-2 



r 




2r + £ + a + e-2 
r-1 




z 



,2r+2^+a+e-2^£ 



Now we apply (2.7) to the inner summation and obtain 



1 - KZ^C{Z^) 



1 



1- y/l- 
2z 




(2.11) 



Summing the expressions (2.9) and (2.11) gives exactly (2.2). □ 

3. The combinatorics of n vicious walkers. In this section we recall Lind- 
strom's classical result [38, Lemma 1] on the enumeration of families of vicious walkers 
(non-intcrsccting paths)^ in directed graphs (directed networks), and apply it to ob- 
tain determinantal formulae for our partition function zj:^\a. — > e; k). 

Let G = (y, E) be a directed acyclic graph with vertices V and directed edges E. 
Furthermore, we are given a function w which assigns a weight w{x) to every vertex or 
edge X. Let us define the weight w{P) of a walk P in the graph by Yle w(e) Ylv 'w{v), 
where the first product is over all edges e of the walk P and the second product is 
over all vertices v of P. We denote the set of all walks in G from to f by V{u — > v), 
and the set of all families (Pi, P2, • • • , -Pn) of walks, where Pi runs from ttj to Vi, 
i = 1, 2, . . . , n, by V{u v), with u = (tti, tt2, . . . , Un) and v = {vi,V2,..., Vn). The 
symbol V~^{u v) stands for the set of all families (Pi, P2, . . . , Pn) in V{u v) 
with the additional property that no two walks share a vertex. We call such families 
of walk(er)s "vicious walkers" or, alternatively, "non-intersecting lattice paths". The 
weight w(P) of a family P = (Pi, P2, . . . , P^) of walks is defined as the product 
nr=i ^(-f'l) of all the weights of the walks in the family. Finally, given a set M with 
weight function w, we write GF{A4;w) for the generating function XlxeA^ w{x). 

We need two further notations before we are able to state the Lindstrom-Gessel- 
Viennot theorem.^ The symbol denotes the symmetric group of order n. Given a 

^Lindstrom used the term "pairwise node disjoint paths". The term "non-intersecting," which is 
most often used nowadays in combinatorial literature, was coined by Gessel and Viennot [24]. 

^By a curious coincidence, Lindstrom's result (the motivation of which was matroid theory!) was 
rediscovered in the 1980s at about the same time in three different communities, not knowing from 
each other at that time: in statistical physics by Fisher [17, Sec. 5.3] in order to apply it to the 
analysis of vicious walkers as a model of wetting and melting, in combinatorial chemistry by John 
and Sachs [30] and Gronau, Just, Schade, Scheffler and Wojciechowski [28] in order to compute 
Pauling's bond order in benzenoid hydrocarbon molecules, and in enumerative combinatorics by 
Gessel and Viennot [24, 25] in order to count tableaux and plane partitions. Since only Gessel and 
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permutation u e /Sn, we write u^- for (■Uo-(i)) '^a{2)j • • • ) '^a{n))- Then 




Most often, this theorem is apphed in the case where the only permutation a for 
which vicious walks exist is the identity permutation, so that the sum on the left- 
hand side reduces to a single term which counts all families (Pi, P2, ■ ■ ■ , Pn) of vicious 
walks, the i-th walk Pi running from Ai to Ei, i = 1, 2, . . . , n. This case occurs for 
example if for any pair of walks (P, Q) with P running from Ua to and Q running 
from Ub to Vc, a < b and c < d, it is true that P and Q must have a common vertex. 
Explicitly, in that case we have 

GF{V+{u^v);w) = det (GF{V{uj^Vi);w)). (3.2) 

l<i,j<n 

This is also the case which we encounter in this paper. ^ 

Proposition 2. Let a — (ai, 02, ... , a^) and e = (ei, 62, . . . , e^) be n-tuples of 

non-negative integers with < ai < a2 < ■ ■ ■ < an, < ei < 62 < ■ ■ ■ < e^, all ai's 
of the same parity, all 's of the same parity, such that Ui + Ci = t (mod 2) for all 

i. As before, let zj:^\a e; k) denote the partition function for families of n vicious 
walkers, the i-th starting at (0, a^) and ending at {t, e^), none of them running below 
the X-axis, where the weight of a walker configuration P is defined as n'^^^^ with c(P) 
denoting the number of contacts of the walkers with the x-axis. Then 

Z[''\si ^ e; k) = det (z^^\aj ef, k)) , (3.3) 

l<i,j<n \ / 

where Z^^\a — > e; k) is given by (2.1). 

Proof. We apply (3.2) for G the graph with vertices the points {x, y) in the plane 
integer lattice with y > and with directed edges {x,y) —>■ {x + l,y — 1), y > 1, and 
(x, y) — (x + 1, y + 1) , y > 0. As the weight w we choose the function which assigns 1 
to every edge and to every vertex (x, y) with y > 0, and which assigns k to every vertex 
on the a;-axis. If we now choose Ui — (0, ai) and Vi = (t, e^) in the Lindstrom-Gessel- 
Viennot theorem (3.2), then the generating function on the left-hand side is exactly 
the partition function zj:"'^ (a — > e; k) for n vicious walkers, whereas the generating 
function which gives the (i,j)-entry of the determinant on the right-hand side is the 



Viennot rediscovered it in its most general form, I propose to call this theorem the "Lindstrom— 
Gessel— Viennot theorem." It must however be mentioned that in fact the same idea appeared even 
earlier in work by Karlin and McGregor [32, 33] in a probabilistic framework, as well as that the 
so-called "Slater determinant" in quantum mechanics (cf. [48] and [49, Ch. 11]) may qualify as an 
"ancestor" of the Lindstrom— Gessel— Viennot determinant. 

^There exist however also several interesting applications of the general form of the Lindstrom— 
Gessel- Viennot theorem in the literature, see [10, 16, 51]. 
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partition function Z^\aj — > e^; k) for a single walker from (0, aj) to (t, e^), the latter 
being given explicitly by Proposition 1. □ 

This result has also been obtained by Brak, Essam and Owczarek in [9, Eqs. (22), 
(23)] (using the Bethe Ansatz method; their derivation is much more complicated, 
but it is on the other hand a method which is more widely applicable that they use). 

In this paper we will primarily analyse the case of watermelon configurations which 
start on the wall. To be precise, this is the case where the starting point Ui is (0, 2z — 2) 
and the end point is (t, y + 2z — 2), z = 1, 2, . . . , n, for some non-negative integer y. 
In our analysis of this particular case we shall in fact use two variants of the above 
formula, which are also corollaries of the Lindstrom-Gessel-Viennot theorem. 

Proposition 3. Let t and y be non-negative integers with t = y (mod 2). As in 
the Introduction, let zj:^\y; k) denote the partition hmction for famihes ofn vicious 
walkers, the i-th starting at (0, 2i — 2) and ending at {t,y -\- 2i — 2), i = 1, 2, ... , n, 
none of them running below the x-axis, where the weight of a walker configuration P 
is defined as n^^^^ with c(P) denoting the number of contacts of the walkers with the 
X-axis. Then 

where 

_^y + 2i + £-!/' t + 2j-i \ ^ 



(3.5) 



Furthermore, let Z^r (k) denote the partition function for families ofn vicious walkers, 
the i-th starting at (0, 2i — 2) and ending at (2r, 2i — 2), i = 1, 2, . . . , n, none of them 
running below the x-axis, with the same weight of walker configurations. Then 

^2?(^) = ^i2n<d.|t (C{r + i+j;K)), (3.6) 

K 0<i,j<n—l 

where 



Proof. We start by proving the first claim. Given a family of vicious walkers as 
in the first statement of the proposition, we may freely attach 2i — 2 up-steps at the 
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Watermelon configuration with attached walk pieces at the beginning 



Figure 4 



beginning of the i-th walk, i = 1,2, ... ,n, see Figure 4 for the resulting walks if we 
do this with the watermelon configuration in Figure l.b. 

It is obvious that the number of families of vicious walkers with starting points 
u[ — {—2i + 2, 0) (instead of Ui = (0, 2i — 2)) and end points Vi = {t, y + 2i — 2) , i = 
1,2, ... ,n, each walker not running below the x-axis (see Figure 4), is exactly the same 
as the number of families of vicious walkers in the first statement of the proposition 
(compare Figure 1). Clearly, by attaching these walk pieces at the beginning, we 
introduced n — 1 further contacts of the first walk with the a;-axis (namely in the 
points u'^, . . . , U2). 

If we now apply (3.2) with Ui replaced by u'-, then we obtain that the partition 
function for the latter families of walkers is given by 

det ( GF(P((-2j + 2, 0) ^ {t, y + 2i- 2)); w)), 

l<t,j<n 

where w is our contact weight. By definition, the generating function 

GF(P((-2j + 2, 0) ^ {t, y + 2i- 2)); w) 

is equal to zj:^2j-2i^ — > y + 2z — 2; k), which by Proposition 1 is exactly the same as 
Bi-ij-i{t, y; k). Since our operation of adding walk pieces introduced n — 1 additional 
contacts with the x-axis, we must divide the above determinant by k""^, and we 
obtain the final result (3.4) after replacing i hy i + 1 and j by j ' + 1. 

The proof of formula (3.6) is completely analogous. Here, one has to also attach 
2z — 2 down-steps at the end of the z-th walk. We leave the details to the reader. 
□ 
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4. Exact results for the partition function for watermelons of devia- 
tion 0. In this section we address the exact evaluation of the partition function of 
watermelon configurations which start and end on the wall. To be precise, we con- 
sider the partition function Z^\k) = Xlp^'^''^''' where the sum is over all families 
P = [Pi, P2, . . . , Pn) of n vicious walkers never running below the x-axis, the i-th 
walker P^ starting at (0, 2i — 2) and ending at (2r, 2z — 2), i = 1, 2, . . . , n, and with 
c(P) denoting the total number of contacts of the walkers with the wall. It should 
be noted that Z^^^(k) = z£^((0, 2, 4, . . . , 2n - 2) {0,2,4, ... ,2n - 2); k) in the 
notation of the Introduction and of Proposition 2. 

The theorem below provides an expression for the partition function [k] in 
terms of a single sum, by finding a product formula for the number of the above 
watermelon configurations with exactly I contacts with the wall. This result was first 
obtained by Brak and Essam [6, Theorem 6] by clever determinant manipulations, 
which allowed them to find a recurrence for the partition function. Our proof is 
completely different. It demonstrates that, once we have a guess for the formula 
(which can be found completely automatically using the Mathematica program Rate^, 
respectively its Maple equivalent GUESS^), it is already proved, as the validity of the 
formula can be established in a completely automatic fashion, by making use of what 
is known under the name of "Dodgson's condensation method" and the Gosper- 
Zeilberger algorithm. A more general result will be obtained later in Theorem 8, 
however using completely different (non-automatic) methods. 

Theorem 4. The partition function (k) is equal to 

jr - 1)! nr=o (2^ + 1)! nr=o'(2r + 2z)! g - £ - 2\ /£ + 2n - 1\ ^^^^^ 



e=o 



In particular, the number of walker configurations with exactly i contacts with the 
X-axis is equal to 

(r - 1)! nL"o (2^ + 1)! nL"o'(2r + 2z)! f2r - A /£ + 2n - 3 



r-l \ 1-2 



Proof. "Dodgson's condensation method" (cf. [36, Sec. 2.3]) is based on the fol- 
lowing determinant identity due to Desnanot and Jacobi. Let A be an n x n matrix. 
Denote the submatrix of A in which rows ii, 12, . ■ ■ ,ik and columns ji, J25 ■ ■ ■ ^jk are 
omitted by A-'-'^'p''"'-'' . Then there holds 

det A ■ det Alll = det A] ■ det A"^ - det ■ det A^. (4.2) 

If n = 2, then det A\'^, the determinant of a x matrix, has to be interpreted as 1. 



^written by the author; available from tittp : //www . mat . univie . ac . at /~kratt 



^written by Francois Beraud and Bruno Gauthier; available from 



littp : //www-igm .univ-mlv . f r/~gauthier 
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Because of Proposition 3, Z2^\k) is essentially equal to a determinant. We now 
apply the above determinant identity with A equal to the determinant in (3.6). Thus 
we obtain 

Zir\^)zi:^f{^) = Zi-l\^)zir'\^) - {Zi-i\^)y- (4.3) 

This relation enables us to prove the theorem by induction on n. It is clearly true for 
n = 1, thanks to Proposition 1. Equation (4.3) allows to perform the induction step. 

We only have to check that (4.3) is true with Z^^\k) as asserted in the statement of 
the theorem. 

Let first n = 2. Then the term Z'2^^4(k) appears in (4.3), which by virtue of (3.6) is 
times the determinant of a x matrix. Since, according to our earlier convention, 

the latter determinant should be interpreted as 1, we should let Z^^I\k) := k^. Thus, 

for n = 2 Equation (4.3) becomes 

r(r + l)!(r + 2)!^ V r-1 

_ /y^ £j-^/2r-£ + 2\ ^+2^ fy^i + l f2r-i-2 



- E 



. f + 1 



By comparison of coefficients of k^^'^ on both sides, this is seen to be equivalent to 
verifying that 



(e+l)(e + 2)(e + 3) 

e 

= E 



2r - e - 2 
r — 1 

^ + 1 /2r - £ + 2\ e - £ + 1 /2r - e + £ - 2 



^=0 



r-\-2\ r+1 J r \ r — 1 

e+1 /2r-e\e-e + l /2r-e + e 



r + 1\ r J r + 1 

This identity is easily proved once one observes that the summand of the sum on the 
right-hand side is equal to G(e,£ + 1) — G(e, ^), where 

^ (2r-£+l)! (2r-e + e-2)\ 

^ ' ' (r + 1)! (r + 2)! (r - £ + 1)! (r - e + £ - 1)! 

X (-6 - 8e - 2e^ + U + Gel + e^i -f + 2ef + e^f - - 2e^^ 

+ e^- 6er - 2eV - Sir - Aelr - e^lr + Zfr + e£V + Gr'^ + 2er'^), 
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since then the sum on the right-hand side is equal to G{e,e + 1) — G(e, 0), which 
can be seen to be equal to the left-hand side upon some simplification. (Clearly, 
Gosper's algorithm [26], [27, §5.7], [43, §11.5] was used to find G{e,£). The particular 
implementation that we used is the Mathematica implementation by Paule and Schorn 
[42].) 

Now let n > 2. In this case, substitution of the claimed expression for Zg"^ (k) in 
(4.3) and comparison of coefficients of k'^'^'^ on both sides yields, after some manipu- 
lation, that we have to verify the summation 



E 



(n - 1) (2n - 1) (r + 1) /2r - £ - 2\ /2n + £ - 1 
(2n + r-2)(2r + l) V r-1 )\ £ 

'2r -e + e + 2\ f2n + e - £ - 5 

(r + 1) (n + r - 1) (2n + 2r - 3) /2r - £ + 2\ /2n + £ - 3 
(2n + r - 2) (2r + 1) V ^ + 1 /V 
'2r - e + £ - 2\ /2n + e - £ - 3 
r-1 )\ e-i 



2r 



2n + £ - 3\ /2r - e + A /2n + e - £ - 3' 
i )\ r )\ e-i 



= 0. (4.4) 



The most straight-forward way to do this is to feed the sum into the Gosper-Zeilberger 
algorithm [43, 54, 55] (again we used the Mathematica implementation by Paule and 
Schorn [42]). If S{n) denotes the sum on the left-hand side of (4.4), then the output 
of the algorithm is the recurrence 

(n - 1) (2n - 1) (e + 4n - 5) (2n + r - 1) (2n + r) S{n + 1) 

- (e + 4n - 1) (e + 2n - r - 2) (e + 2n - r - 1) (n + r) (2n + 2r - 1) S{n) = 0. 

(4.5) 

Now, it is straight-forward to check that S{1) = 0, thus proving (4.4) for n= 1. The 
validity of (4.4) for arbitrary n then follows upon induction on n from the recurrence 
(4.5). 

Thus, the proof of the theorem is complete. □ 

If we reverse the order of summation in the sum in (4.1) (i.e., if we replace £ by 
r — 1 — £) and then rewrite it using the standard hypergeometric notation 



p-Tq 



ai , . . . , a 
bi,...,bq 



'P . 



= E 



{ai)m ■ ■ ■ jap) 
ml {bi)m ■ ■ ■ {bg) 



where the Pochhamnier symbol {a)m is defined by (a), 
m > 1, {a)o '■= 1, then we obtain 



a{a + 1) ■ ■ ■ {a + m — 1), 



{r-l)V 



r, 1 — r 1 
2 - 2n - r ' K 
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or, more precisely (the above 2-Pi-series is actually ill-defined because of the denomi- 
nator parameter 2 — 2n — r which is a negative integer), 



lim K 



l+r 



(2n) 



r-l 



(r- 1)! 



r, 1 
2-2n- 



■ r 1 



(4.6) 



In view of the fact that there are numerous 2-f^i-transformation formulae, we can 
obtain numerous equivalent expressions for this sum and, hence, for (4.1) itself. 

A particular such expression has been given by Owczarek, Essam and Brak in [41, 
Eq. (4.65)], with proof provided by Brak and Essam in [7], and was the starting point 
for their scaling analysis of the partition function Z2^\k) in [41]. This expression 
provides an expansion of (4.1) as an (infinite) power series in {k — I)/^^. It is however 
only valid for k < 2. In the corollary below, we show how the expression of Owczarek, 
Essam and Brak can be readily derived by applying appropriate 2i^i-transformation 
formulas to (4.6). Moreover, this technique yields also an expression which is valid for 
K > 2. In the case n = 1 the latter expansion was, for example, previously observed 
by Brak and Essam [5, Eq. (25)]. 

Corollary 5. Let k > 2(-\/2 — 1). Then the partition function is equal to 



niro'(2z + l)!(2r + 2z)! 



1 



K 



2n-2 



E 

h=0 



n + h 
h 



1\ (2r + 2n-r 



2h 



(r + n)h (r -|- 2n)h \ k 



K 



1 



+ x(ac>2) 



(k - 2)K2r+2n-i (r + 2n-l)! 



(k-1)^+2"-i (2n- 1)! (r + n- 1)! 

"""^ 'n-l\{r + n-h-l)\{r + 2n-h-2)\ 
h 



E 



(2r + 2n-2h-2)\ 



K 



(4.7) 



where x(^) = 1 if A is true and x(^) = otherwise. 

Proof. First let k < 2. We write the sum in (4.1) in the form (4.6) and apply 
the transformation formula (cf. [50, (1.8.10), terminating form]) 



a. 



-N 



(c-a) 



N 



(c) 



N 



-N,a 
c - 



where N is a non- negative integer. Now the limit e 
we obtain 



1+a-c- N' 

can be safely performed, and 



K, 



^^^{2n + r)r-i 



(1) 



-2-fl 



r-l 



1 



r, r _ ^ _ i. 
2n + r ' K 



for the sum in (4.1). Next we apply the quadratic transformation formula (see [44, 
(3.31), reversed]) 



2-^1 



a, 1 



= {l-z) 



c-1 



2-ri 



2 2' 2^+ 2 2.4^(1_^) 



(4.8) 
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which holds when Hez < 1/2, as well as if 2; = 1/2, and substitute the result back 
in (4.1). After some manipulation, we obtain the expression in the first line of (4.7), 
which, taking into account the condition under which (4.8) is true and the latter series 
converges, holds for 2(^2-1) <k<2. 

If on the other hand k > 2, then we proceed differently. We start with the expres- 
sion 



lim K, 



l+r 



(2n) 



r-l 



(r- 1)! 



r — e,l — r + e ^ 1 
2-2n-r + £''^ 



= lim 



{2n)r- 



e^o (r-l)! ^ (2-2n-r + £)fcA;! 



(4.9) 



which is clearly equivalent to (4.6). The reader should however note the differences 
to (4.6) in placing the e. We split the sum over k into the three ranges < /c < r — 1, 
r<k<r + 2n — 2, and r + 2n — 1 < k, and then we perform the limit £ — > 0. 
Thereby each summand in the range r<k<r + 2n — 2 vanishes identically. Thus, 
upon writing the sum over the range r + 2n — 1 < A; in hypergeometric notation, the 
expression (4.9) becomes 



K 



l+r 



(2n 
(r 



r-l 



tjr-1 (r)fc (1 - r)fc 

-l)!^(2-2n-r)fc/c 



(r) 



r+2n-l 



^2n-2 (^ + 2n-l)! 



2n, 2r + 2n-ll' 
2n + r ' K 

(4.10) 

The first sum in this expression is exactly the sum in (4.1) (which is most easily seen 
from the equivalent expression (4.6)). Therefore, if we equate (4.9) and (4.10), we see 
that the sum in (4.1) is equal to 



(r) 



r+2n-l 



^2n-2 (r+2n- 1)! 



2-^1 



2n, 2r + 2n - 1 J- 
2n + r 'k 

+ lim K 



l+r 



(2n) 



r-l 



(r-l)! 



r — £, 1 — r + £_ 1 
2-2n-r + £'K 



To the first 2-^1-series we apply the quadratic transformation formula (cf. [4, Ex. 
4.(iii), p. 97, reversed]) 



iFi 



a, h 

1_ I a I 6 ) ^ 
2 2 2 



— 2-^1 



2 2 2 



which is valid for Rcz < 1/2, while to the second 2-^1 -series we apply the quadratic 
transformation formula (4.8), and subsequently the transformation (cf. [50, (1.3.15)]) 



2^1 



a, h 



(1 - 2)^-"-^Fi 



c — a, c — 6 



If the resulting expression is simplified and then substituted back in (4.1), we obtain 
eventually the expression (4.7) for k > 2. □ 
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Remark. (1) Other alternative expression for (4.1) will be obtained later (in 
greater generality) in Theorems 9 and 11. 

(2) The sum in (4.1) does not simplify in general. However, there are two cases 
where it does, namely if k, — 1 and if k, — 2. We state these results (again in greater 
generality) in Corollary 10. 

5. Vicious walkers, semistandard tableaux, and jeu de taquin. The pur- 
pose of this section is to show that the number of watermelon configurations which 
start on the wall, end at some deviation y, and have a given number of contacts with 
the X-axis is equinumerous to another family of vicious walkers which do not run 
below the but have no restriction on the number of contacts with the 

The upshot of this is that the latter number of watermelons can now be written in 
form of a single determinant by applying the Lindstrom-Gessel-Viennot determinant 
(3.2). The result has been originally found by Brak and Essam [6, Cor. 1]. They 
derived it by manipulating (a variant of) the Lindstrom-Gessel-Viennot determinant 
(3.4). Since the result itself is purely combinatorial, they posed the problem of finding 
a bijective proof. We resolve this problem here. The bijective proof that we shall give 
below is based on the fact that vicious walker configurations are in bijection with 
semistandard tableaux (see [29, 37]) and on a modified jeu dc taquin for semistan- 
dard tableaux that occurred earlier in a bijective proof [35] of Stanley's hook-content 
formula for the number of semistandard tableaux of a given shape with bounded 
entries. 

Proposition 6. The number of families of n vicious walkers, the i-th starting at 
(0, 2i — 2) and ending at {t,y + 2i — 2), i — 1,2, ... ,n, none of them running below 
the X-axis, and where the first walk has £ + 1 contacts with the x-axis {including its 
starting point (0,0)) is the same as the number of families of n vicious walkers, the 
i-th starting at (0, 2i — 2) and ending at {t,y + 2i — 2), i = 1, 2, . . . , n — 1, the n-th 
walker running from (0, 2n — 2) to {t — £ — l,y + 2n + £ — 3), none of them running 
below the 

Remark. The statement of the proposition deviates slightly from the one in [6, 
Cor. 1]. The equivalence of both statements can be seen by observing that in the 
formulation of Brak and Essam there are forced walk portions at the beginning and 
at the end of the walks, which can hence be omitted. There is also a typo there: the 
definition of v'^ must be {t — m,y + 2{n — 1) + m,). 

Proof of Proposition. It is well-known that configurations of vicious walkers 
are in bijection with so-called semistandard tableaux. By definition, a semistandard 
tableau is an array of integers in which entries are weakly increasing along rows and 
strictly increasing along columns. Figure 5 shows three examples. 

To explain the bijection, let us first consider the watermelon configurations in the 
statement of the proposition (this is the first set of families of vicious walkers in the 
statement of the proposition). An example for n = 4, t = 12, y — 2, and £ = 3 is 
shown in the Figure 6. a. Now label down-steps by the x-coordinate of their starting 
point, so that a step from (a, b) to {a + 1, b — 1) is labelled by a, see Figure 6. Then, out 
of the labels of the z-th walk, form the i-th. column of the corresponding tableau. If 
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b. 

Semistandard tableaux 
Figure 5 
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c. 



this procedure is applied to the family of vicious walkers in Figure 6. a then the tableau 
in Figure 5. a is obtained. It is not difficult to see that the resulting array of numbers 
is always a semistandard tableau. In fact, the entries are trivially strictly increasing 
along columns, and they are weakly increasing along rows because the walks do not 
touch each other. Moreover, the restriction that no walker may run below the x-axis 
translates into the condition that the entries in the i-th. row of the corresponding 
tableau are at least 2i — 1. Finally, any contact of the first walker with the a;-axis 
(except the one in (0, 0)) produces an entry in the first column which is exactly equal 
to its lower bound. (In the tableau in the Figure 5. a these are the entries 1, 7 and 9.) 

Thus, the first set of families of vicious walkers in the statement of the proposition 
is in bijection with semistandard tableaux with n columns of length (t — y)/2, with 
the entries in row i at least 2i — 1 and at most t — 1, i = 1, 2, . . . , (t — j/)/2, and with 
exactly £ entries in the first column which are equal to their respective lower bound. 

If we apply this same translation to the second set of families of vicious walkers 
in the statement of the proposition, then we see that these are in bijection with 
semistandard tableaux with n columns, the first n — 1 of them of length {t — y)/2, 
the last of length {t — y — 2i)/2, with the entries in row i being at least 2i — 1 and 
at most t — 1, i = 1,2, . . . ,{t — y)/2, and with the entries in the last column being 
at most t — £ — 2. Figure 5.c shows the result when this translation is applied to the 
family of vicious walkers in Figure 6.b. 

It is now our task to construct a bijection between the latter two sets of tableaux. 
We start with a tableau from the first set (of which an example is shown in Figure 5. a). 
The construction makes use of a modified form of jeu de taquin from [35, (2.2)-(2.4)]. 
We start with the bottom-most entry in the first column which is equal to its lower 
bound. (In Figure 5. a this is the entry 9 in the first column.) Let us call it the 
"special entry," and denote it by s. Now we move s to the right/bottom by the 
following procedure, denoted by (JT): 

(JT) Compare the special entry s with its right neighbour, x say (if there is no 
right neighbour, then, by convention, we set x = oo), and its bottom neighbour, y 
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b. 

Corresponding vicious walkers 



Figure 6 



say (if there is no bottom neighbour, then, by convention, we set y — oo), see (5.1). 
If there is no right or bottom neighbour, then stop. 
If not, then we have the following situation. 



s 


X 


y 





(5.1) 



where at least on of x and y is not oo. If a; > y — 1 then do the move 





X 


s 





li x+1 < y then do the move 



X + 1 


s 


y 





(5.2) 



(5.3) 



Repeat (JT). 

For example, if we apply this algorithm to the tableau in Figure 5. a, the special 
entry being the 9 in the first column, then we obtain the tableau in Figure 7. (In this 
case it is only movements of type (5.3) that are applied.) 

Clearly, after this algorithm has terminated, the special entry will have ended up in 
the bottom-right corner (as it does in this example). We now remove the special entry 
(the 9 in the bottom-right corner in our example). 
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Figure 7 



Next, the same algorithm is apphed to the (now) bottom- most entry in the first 
column which is equal to its lower bound. (In Figure 7 this is the entry 7 in the first 
column.) Since successive jeu de taquin paths cannot cross each other (see [35, Fig. 9] 
and the accompanying explanations for a detailed argument), this entry must again 
end up in the bottom of the last column. Subsequently, we remove it. This procedure 
is repeated until there is no entry in the first column which is equal to its lower bound. 
(The final result when this procedure is applied to the tableau in Figure 5. a is the 
tableau in Figure 5.b.) 

Now, what do we obtain in the end? It is easy to see that after each movement of 
type (5.2) or (5.3) the entries are weakly increasing along rows and strictly increasing 
along columns, if we ignore the special entry. Thus, we will have obtained a semistan- 
dard tableau with n columns, out of which all of them have length {t — y)/2 except 
for the last, which will be by £ shorter. Furthermore, all entries will be at most t (this 
is due to the fact that at most once 1 is added to an entry by performing a movement 
of type (5.3)), and all entries in the z-th row will now be strictly larger than 2z — 1. 
Finally, since originally the entry in row i in the last column is at most (t+y + 2z — 2)/2 
due to strict increase of entries along columns, and since whenever an entry (other 
than the special entry) is moved upward (which is only possible through a movement 
of type (5.2)) it shrinks by 1, the entries in the last column will in the end be at most 
t — I — 1. Clearly, if we subtract 1 from all the entries of such a tableau, we obtain a 
tableau from the second set of tableaux. (In our running example, after subtracting 1 
from all the entries in the tableau in Figure 5.b we obtain the tableau in Figure 5.c.) 

We claim that this is a bijection between the sets of tableaux under consideration. 
In order to see this, we have to construct the inverse mapping. In fact, every single 
step can be inverted. We would start with a tableau from the second set, add 1 to 
all the entries, place an entry s (whose value will be determined later) immediately 
below the bottom-most entry in the last column of the tableau, and then apply the 
following "inverse" jeu de taquin, denoted by (JT*): 

(JT*) Call s the special entry. 

If the special entry s is located in the first row and column, then stop. 
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If not, then we have the following situation, 





y 


X 


s 



where one of x or y could also be absent. (If y is actually not there, then by convention 
we set y = — oo.) We denote the number of the row in which s is located by i. 

If X is actually not there and y + I < 2z — 1, then stop. Otherwise, and also if 
X < y + 1, then do the move 





s 


X 


y + 1 



If X — 1 > y then do the move 





y 


s 


x-l 



Repeat (JT*). 

Suppose that the special entry ended up in row i. Then we now set s = 2z — 1. 

Subsequently, this procedure is repeated until the complete {{t — y)/2) x n rectangle 
is filled. We leave it to the reader to check that this yields indeed the desired inverse 
mapping. □ 

Remark. Brak and Essam have also a variant of Proposition 6 in [6, Cor. 2]. It 
states that the number of families of n vicious walkers, the z-th starting at (0, 2i — 2) 
and ending at {t,y + 2i — 2), i = 1,2, ... ,n, none of them running below the a;-axis, 
and where £ of the contacts of the first walk with the x-axis (excluding (0,0)) are 
marked, is the same as the number of families of n vicious walkers, the i-th starting 
at (0, 2z — 2) and ending at {t, y + 2i — 2), i = 1,2, . . . , n — 1, the n-th walk running from 
(0, 2n — 2) to {t, y + 2n + 2£ — 2), none of them running below the a;-axis. Again, they 
prove it by determinant manipulations and pose the problem of finding a bijective 
proof. Such a bijective proof can be provided in the same manner as the one that 
we found for Proposition 6. First, one translates both sets of walk families into sets 
of semistandard tableaux in the same manner as before. The first set translates into 
semistandard tableaux with n columns of length {t — y)/2, with the entries in row i at 
least 2i — l and at most t — 1, i= 1,2, . . . , {t — y)/2, and with exactly £ entries in the 
first column which are equal to its lower bound and which are marked. The second 
set translates into semistandard tableaux with n columns, the first n — 1 of them of 
length (t — y)/2, the last of length {t — y — 2£)/2, with the entries in row i being at 
least 2i — 1 and at most t — 1, i = 1,2, . . . , {t — y) /2. The bijection of the proof of 
Proposition 6 then becomes a bijection between the latter two sets of tableaux when 
one replaces the modified jeu de taquin (JT) and (JT*) by ordinary ]e\x de taquin (as 
for example explained in [47]; i.e., nothing is subtracted from or added to x and y in 
(5.2), (5.3), (5.4), (5.5); furthermore the condition x < y -\- 1 has to be replaced by 
X < y everywhere, and x — 1 > y has to be replaced by x > y everywhere, as well as 
the stopping rule in (JT*) has to be modified to the extent that y + 1 <2i — l has to 
be replaced hy y < 2i — 1). 
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6. Exact formulas for the partition function for watermelons of arbi- 
trary deviation. This section is devoted to the derivation of exact formulas for 
the partition function zj:^\y; k) for families of n vicious walkers, the i-th starting at 
(0, 2i — 2) and ending at {t,y + 2i — 2), i — 1,2, ... ,n, none of them running below 
the a;-axis, as defined in (1.1). The main result is Theorem 8, which expresses the 
partition function in terms of a double sum. It extends Theorem 4 to watermelon 
configuration of an arbitrary (but fixed) deviation. By applying hypergeometric sum- 
mation and transformation formulas to this double sum, we are also able to provide 
alternative expressions for z[^\y; k), special cases of which having been given earlier 
by Brak and Essam (private communication). 

Our starting point is Proposition 6. We apply the Lindstrom-Gessel-Viennot for- 
mula (3.2) to the second set of vicious walkers in the proposition in order to obtain a 
determinant for the number of watermelons of deviation y and exactly £ + 1 contacts 
with the The most useful determinant is obtained if we first prepend 2i — 2 

up-steps to the i-th walk (analogously to what we did in the proof of Proposition 3) , 
so that we obtain families of n vicious walkers, the z-th starting at {—2i + 2, 0) and 
ending dX {t,y -\-2i — 2), i = 1,2, ... ,n — l, the n-th walker running from (— 2n -|- 2, 0) 
to (t — £ — 1, y + 2?i + £ — 3), none of them running below the x-axis. Again, it is clear 
that the prepended walk portions are forced, and, hence, the number of the latter 
families of walkers are the same as the number of families in the second set of walkers 
in Proposition 6. If we apply the Lindstrom-Gessel-Viennot formula to the modified 
families of walkers, then we obtain 

/ y+2i-l ( t+2j-2 X 1 <,•<„- 1 \ 

det f ^+^+^-^ l<i<n 1\ 

l<iT<n ( t+2j-£-3 ^ ■ I ' l^.i; 

where we used (2.3) to compute the number of walks with given starting and end 
point, which do not pass below the x-axis. 

As it turns out, this determinant can be simplified to a single sum. We state the 
corresponding result in the lemma below. 

Lemma 7. Let t and y he non-negative integers with t = y (mod 2). Then the 
number of famihes of n vicious walkers, the i-th starting at (0, 2i — 2) and ending 
at {t,y + 2i — 2), i — 1,2, ... ,n, none of them running below the x-axis, and where 
the Rrst walk has £+1 contacts with the x-axis {including its starting point (0, 0)) is 
equal to 

-izi {t + 2iy.i\ -^ (y + n + z-iy. 
to C-^ + n + i-iy.C-i^ + ty. l{ {y + 2iy 

V\ l^n-/c-l iy + 2n + i-2)it + 2k-i-iyC-i^ + ky 

' {t + 2ky.{^ + k-n-£+iyk\{n-k-iy.' ^ ' ' 

Proof. We expand the determinant (6.1) along the last row. Thus we obtain 



E(-i) 



,+fey + 2n + £- 2 ( t + 2k-£-2> 
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where is the minor of (6.1) where the last row and the k-th column has been 
deleted. This minor is easily evaluated by means of [36, Theorem 26, (3.13)] with 
A = 1 — y, B — 2, and Li = i + x{i > ^) + — 1, where x(A) = 1 if A is true and 
x(A) = otherwise. After some simplification, and after replacing k hy k + 1, one 
obtains (6.2). □ 

The preceding lemma provides us with a workable expression for the number of 
watermelon configurations with a given number of contacts with the x-axis. We will 
use it now to derive several expressions for the partition function Zf^\y; k). The first 
of those will be the one which is most suited for studying the asymptotic behaviour 
of the partition function, and for the mean number of contacts, which we will do in 
Sections 7 and 8, respectively. 

Theorem 8. Let t and y be non-negative integers with t = y (mod 2). The 

partition function zj:^\y; k) for faniihes of n vicious walkers, the i-th starting at 
(0, 2i — 2) and ending at {t,y + 2i — 2), i — 1,2, ... ,n, none of them running below 
the X-axis, where the weight of a walker configuration P is defined as k,'^^^^ with c(P) 
denoting the number of contacts of the walkers with the x-axis, is given by 

C-^-iy. H {t + 2i)\i\{y + n + t-l)\ 



(l+I + 2n-2)! f-J^ {t^+n + i-iy.C-^ + iy.{y + 2iy. 

{t-y)/2n-l 



X 



E E(-i)' 



n-l\/y + 2k-2\/ t-l-l 

. k )[ 2k )y-^-£-k 
e=o k=0 ^ ^ ^ / \ 2 

^ + £ + 2n-2\ (^),(2n-2fc-l)!(2/c)! 



2n-2A;-l ; (i^+^-fc)^ 



(6.3) 



Remark. In the case n — 1 the product in (6.3) is empty and has to be interpreted 
as 1. Moreover, because of the binomial coefficient ("^^^) in the summand, the sum 
over k collapses to just one term, the term for k = 0. 

Proof. Clearly, the partition function zl.^\y; k) is equal to 



(i-y)/2 



(6.4) 



£=0 



where N^'^\y;£) is the expression in (6.2). We write the sum in (6.2) in hypergeo- 
metric notation, thus obtaining 



( ^^n+i (?/ + 2n + £-2)(t-£-l)!(^)! 



{n-iytiC-^ -n-e + iy 

X4F3 



^2 2' 2^2'2 2'2' 



2' 
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To the 4F3-series we apply one of the balanced (which, by definition, means that the 
sum of the denominator parameters exceeds the sum of the numerator parameters by 
1) 4F3 transformation formulas due to Whipple (cf. [23, Appendix (III. 16), q t 1]) 



4-^3 



a, b, c, —N 
e,/,l + a + 6 + c-e-/-7V' 

{a)N {-a -b + e + f)N {-a - c + e + f)N 



X 4Fs 



(e)Ar {f)N {-a -b-c + e + /)jv 

-N, —a + e, —a + /, —a — b - 
a 



c+e + f 



b + e + f,-a-c + e + f,l-a-N 
where AT is a non-negative integer. Thus we obtain 

{y + 2n + e-2) {n)t-n-i (t + 2n - * _^ {y + £)2n-2 



(6.5) 



X 4F3 



i_r) ^ y — 1 f — t -I- y 

J- /(., 2' 2'2 



1 4_ 1 4_ M 1 

2 2 2' 2 



+ f,l-n- 



2^2 



y 1 



for the sum in (6.2). If we substitute this back in (6.4), then we obtain the claimed 
expression after some manipulation. □ 

Remark. This theorem generalizes Theorem 4 to watermelon configurations of 
arbitrary (fixed) deviation. Indeed, if we set y = in Theorem 8, then, because of the 
binomial coefficient (^^2fc~^) ^ippearing in the summand in (6.3), the only summands 
which are nonzero in (6.3) are the ones with A; = 0. 

As corollaries of Lemma 7 (which via (6.5) is equivalent to Theorem 8), we are 
able to derive remarkable simple alternative expressions for the partition function 

zj:^\y; k) in terms of single sums. (These expressions however seem to be less suited 
for asymptotic considerations.) We start with an expansion of the partition function 
around k, = 1. 

Theorem 9. With the same assumptions as in Theorem 8, the partition function 
z[^\y] k) is equal to 



n-l 



' n-2 



to J\lUy + 2i)lC-±i^ + n + i)lC-i^+t + l 

{y + 2n + 2h-l){y + 2n + h-2y. 



(t-y)/2 



n + h- 1 

n-l )\y + n + h-iy.C-^ + 2n + h-iy.C-^-hy. 



(6.6) 



Proof. Again, by Lemma 7, the partition function Zf'' [y; k) is equal to (6.4), 
with N^^\y]€) the expression in (6.2). We write k^'^^ = k^^^q ih)^'^ ~ ■'■)'^' 
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substitute this in (6.4). Thus we obtain 



r(n) ^.^ ^n-k-i (y + 2n + £ - 2) (t + 2k - £ - 1)\ C-^ + k)\ 



=0 fc=0 



{t + 2k)\C-^ + k-n-e+l)\k\{n-k-l)\ 



h=0 ^ ^ 



where n^"^-* {y) is the product in the first fine of (6.2). Now we interchange summations, 
so that the sum over i becomes the inner-most sum and the sum over h becomes the 
outer- most sum, and we convert the sum over i to hypergeometric notation. This 
yields 

(t-y)/2„-l 



^r\y)j^, E EM)"-'- 

^ ' h=0 k=0 ^ ' 

(t + 2/c-/i-l)!(^ + /c)! 



(t + 2/c)!(^ + A;-n-/i+l)! 
(y + 2n - 3) 2^1 



/i + l,-| + |- /c + n + /i-l ' 
-t-2A; + /i + l ' 



+ (/i+ 1)2^1 



-t - 2/c + /i + 1 ' 



Both 2-^i-series can be evaluated by means of the Chu-Vandermonde summation 
formula (cf. [50, (1.7.7); Appendix (III.4)]) 



2-^1 



a,-iV 
c ' 



(c — a) 



N 



N 



(6.7) 



where is a non-negative integer. The result is simplified, and the sum over k is 
written in hypergeometric notation. This gives the expression 

nW,,) " (y + 2n + 2/.-l)(ai + „-i)!(ly)! 



3-^2 



l + f + n, |-f + l,l-n 
| + |+n + /i + l,|-f-n-/i + 2' 



The 3F2-series is balanced and can therefore be evaluated by means of the Pfaff- 
Saalschiitz summation formula (cf. [50, (2.3.1.3); Appendix (III. 2)]) 



3^2 



a, b, -N 
c,l + a + b-c- N 



(c - a)jv (c - 6)jv 
(c)iv(c-a-6)7v' 
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where N is a non- negative integer. Subsequent simplification then leads to the claimed 
expression (6.6). □ 

Remarks. (1) The coefficient of k,{k, — 1)'* in (6.6) has a combinatorial interpre- 
tation in terms of watermelon configurations, the i-th walker of the configuration 
running from (0, 2i — 2) to {t,y + 2i — 2), i = 1,2,..., n, but never below the x- 
axis, where contacts with the x-axis other than (0, 0) can be marked or not. It is 
an immediate result of the comparison of (6.6) with the definition of zj:^\y; k). In 
this definition, we replaced k^'^^ by h^^^q (^)('^ — 1)'^- Hence, the coefficient of 
k{k — 1)'^ in (6.6) counts watermelon configurations as above where exactly h contacts 
other than (0,0) are marked. 

(2) It is also possible to prove Theorem 9 in the same (automated) fashion as 
Theorem 4, i.e., by starting with the Lindstrom-Gessel-Viennot determinant (3.4), 
and then arguing by induction, which one would base on the condensation formula 
(4.2), the verification of the binomial summations that one has to establish on the 
way being accomplished by the Gosper-Zeilberger algorithm. 

The sum in (6.6) does not simplify in general. In fact, when written in hypergeo- 
metric terms, it turns out to be a 4F3-series, 



K{y + 2n-l)\ 

74-^3 



(^)!(^ + 2n-l)!(y + n-l)! 



y + 2?1 - 1, I + ?1 + i, 71, + I 



1-K 



|+n-i,y + n,2n+| + | 

(6.8) 

This is a balanced 4F3-series, and there are many identities known for balanced 4F3- 
series, but there is no summation formula available that would apply in full generality. 
However, there are two cases where summation formulas are available, namely if re = 1 
(trivially), and if re = 2. The result corresponding to re = 1, giving the total number 
of the watermelon configurations that we consider in this section (and the subsequent 
ones), has been previously obtained by the author, Guttmann and Viennot in [37, 
Theorem 6]. In fact, in Theorem 6 of [37] it is shown how to derive a more general 
result which is valid for star configurations.^ The result corresponding to re = 2 
has been previously obtained by Ciucu and the author in [11, Theorem 1.4] in an 
equivalent form,^*^ and, in the case y — 0, independently by Owczarek, Essam and 
Brak in [41, Eq. (4.45)]. The proof that we give below is independent from the proofs 
in the papers mentioned. 



®What is shown in [37] is that the result for star configurations is, in equivalent form, part of the 
folklore of combinatorics and representation theory of symplectic groups. A weighted generalization 
of this result, in which each path configuration is assigned a certain q-weight, had been proven 
already in [34, Theorem 7]. 

^''The combinatorial objects that are studied in [11] are rhombus tilings of regions composed of 
equilateral unit triangles. There is a standard bijection between such tilings and families of non- 
intersecting lattice paths, which is, for example, explained in Section 2 of [11], see Figures 2.1, 
2.2, 2.3(a) in that paper. In the particular case addressed by Theorem 1.4 of [11], a 45° rotation 
turns these non-intersecting lattice paths into our vicious walkers. Under these transformations, the 
weight for rhombus tilings defined in [11] becomes our contact weight for vicious walkers, up to a 
mult iplicat i ve const ant . 
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Corollary 10. Let t and y be non-negative integers with t = y (mod 2). Then the 

total number {\)\ 1) of families of n vicious walkers, the i-th starting at (0, 2i — 2) 
and ending at {t,y + 2i — 2), i = 1,2, ... ,n, none of them running below the x-axis, 
is equal to 

{y + n + i)l {t + 2i)l il 



n-l 

n 



The partition function zj:"'\y; 2) for these families of vicious walkers is equal to 

{y + n + i-iy.{t + 2i)\i\ 



(6.9) 



n-l 



L-J {y + 2i)! C-^ + n + i-iy.C-^ + iy.' 



(6.10) 



Proof. The claim (6.9) is immediately obvious once we set k = 1 in (6.6). 

For proving the second claim, we set k = 2 in (6.6) and write the sum in (6.6) in 
hypergeometric notation (equivalently, we set k = 2 in (6.8)). Thus we see that this 
sum is equal to 



2(y + 2n-l)! 



(i^)!(i+M + 2n-l)!(y + n-l)! 



4-^3 



y + 2n- l,f +n + |,n,-| + f 
| + n-|,y + n,2n+| + | 



-1 



The 4F3-series can be evaluated by means of the summation formula (of. [50, (2.3.4.6); 
Appendix (III. 10)]) 



4-^3 



o 
2 ' 



a, 1 + |,6, c _^ 
l + o — 6, l + o — c' 



r(l + a-6)r(l + a-c) 
T{l + a)T{l + a-b-c)' 



Substitution of the result in (6.6) then leads to the claimed expression after some 
simplification. □ 

Our next result provides an alternative expression for zj:^\y; k) in form of an in- 
finite series, which resembles the one in Corollary 5 for the y = case. Again, the 
corresponding formula is only valid for k < 2. As we are going to show, it follows 
readily from the formula in Theorem 9 by applying another well-known hypergeo- 
metric transformation formula. Moreover, by the same techniques we can also obtain 
a similar expression which is valid for k > 2. We presume that, in analogy to the 
analysis in [41], these expressions will be the right starting point in order to carry out 
a scaling analysis of z[^\y; k) near the critical point k = 2. 

Theorem 11. Under the assumptions of Theorem 8, the partition function 
zj:^\y; k) is equal to 



n 

i=0 



(y + n + i-iy. 



(y + 2z)! 




{t + 2iyi\ 



X 



(2-^) 



it C-^ + iyC-^ + n + i-iy^ 

{y + 2n + 2h - ly C-^ + n + h - ly 



K^^'+y-^ f^^^h\{y + n + h-iyC-^ + 2n + h-iy 



K 



(6.11) 
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if 2(-\/2 — 1) < K < 2, and it is equal to 



X 



1 {y + 2n + 2h-l)\{'-^ + n + h-l)\ ( k - l^" 

1 



K'^^+y-^ ^ + 2n + /i - 1)! (y + n + /i - 1)! /i! V «^ 
^ (k- ^ /i! (n-/i-l)!(t + 2n- 2/1-2)! 



«2 



(6.12) 



if K > 2. 



Proof. We work with the expression (6.6) for the partition function z[^\y;n). 
In hypergeometric notation, the sum in (6.6) is equal to (6.8). 

Now let first k < 2. Then to the 4F3-series in (6.8) we apply the quadratic trans- 
formation formula (cf. [4, Ex. 6, p. 97, reversed]) 



4-^3 



a, I + 1, 6, c 
f,l + a-6, 1 + a 



; z 



(1 + ^) 
(1-^)'^+^' 



F2 



^ + ^,l + ^,l + a-6-c, 4z 
l + a-b,l + a-c ' {1- zf 



(6.13) 



which is valid provided l-z] < 1 and 14^/(1 — 2)^ 
K < 2, the sum in (6.6) is equal to 



< 1. Thus, as long as 2(^/2 - 1) < 



(2 



{y + 2n-l)l 



^y+2n-i ^t_yy_ (i±y + 2n - 1)! {y + n- 1)! 



X 3F2 



f+n, |+n+i,n+| + |_ 4{k- 
2n+^ + ^,n + y ' 



This expression is now substituted back in (6.6). Some simplification then yields the 
expression (6.11). 

Next we address the case k, > 2. Clearly, we cannot follow the same line of argument 
because we would not be allowed to apply (6.13) because of |1 — «;| > 1. Instead, we 
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start with the expression 

k{k- if-y^/^ (n)(i_^)/2 {2n + t- l)iy-t)/2 



hm 

£1^0 
£2^0 



X 4-F3 



1 -2n 



2 -n 



t ^ - 
- ^ 1 

2' ^ 



n 



^ 1 

2' 



n 



- - ^ + £9 

9 9 ^ '-^7 



n 



| + |-£2,2-2n-|-| 



£1 



hm 

£1^0 
£2^0 



k,{k — 1) 



(t-2/)/2 



(n)(t-y)/2 (2n + t- l)(y_t)/2 



1)! 



X 



E 

A;=0 



(1 - 2n - t + 2k) (1 - 2n - t)fe 
(1 - 2n - t) k\ 



(l-n-|-|+£2)fc(-| + 



(l-n-| + |-£2)fc(2-2n-|-f-£i)fc Vl-^y " ^^'^"^^ 
We spht the sum over k into the three ranges < A; < + 1 < < + 2n — 2, 

and + 2n — 1 < A;, and then we perform the hmit £1 — > 0. Thereby each summand 
in the range + 1 < A; < + 2n — 2 vanishes identicaUy. Thus, upon writing 



the sum over the range 



t+y 



2n — 1 < /c in hypergeometric notation, the expression 



(6.14) becomes, after also performing the hmit 62 0, 



k{k — 1) 



{n){t-y)/2 {2n + t- l)(y-t)/2 



{t-y)/2 

E 

fc=0 



(1 -2n-t + 2A;) 
(1 - 2n - t) 



(l-2n-t)fe(l-n-|-|)fc(-| + 



K 



A;!(l-n-| + f)fc(2-2n-| 
(¥ + l),+2.-2(^ + 2n-l)^ 



1 



K 



2 "T 2 



1)! 



X4F3 



' ' 2 ^ 2 ' ^ ' 2 2 

2 



1 



1 



K 



(6.15) 



(^_l)?/+2— 1 (£+£ + 2n-l)!(y + n 

■-l + 2n + y, ^ 

-i +71+ |,n + 2/,2n+ I + 2 

As it turns out, the first sum in this expression is exactly the sum in (6.6) (which is 
seen by replacing h by — A; in the sum in (6.6)). Therefore, if we equate (6.14) 
and (6.15), we see that the sum in (6.6) is equal to 

K (y + 2n-l)! 



{k - l)«+2"-i (i=^)! (m + 2n - 1)! (y + n - 1)! 



X 4F3 



-l + 2n + y, I +n + f,n, -| + I 



1 



-| +n+ f ,n + y,2n+ I + I ' 1 - k. 



+ lim 

ei-»0 
£2^0 



AC (« - 1)(* (n)(t_,)/2 (2n + t - l)(,_,)/2 



(^)!(^ + n-l)! 



X 4-F3 



1 -2n 



2 -n 



n 



^ 1 

2' 



n 



I +£2, 



2 2' 2 2 

|,1 - n- I + f - £2,2-2n 



^ + # + £1 



y 

2 ' 2 

^ - ^ - £l 

2 2 
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Now we apply the transformation formula (6.13) to both 4F3-series. Subsequently we 
perform the limit in the second term, and we simplify the resulting expressions. After 
this is substituted back in (6.6), this yields exactly the expression (6.12). □ 

7. Asymptotic formulas for the partition function for watermelons of ar- 
bitrary, but fixed deviation. In this section we embark on the asymptotic analysis 
of the partition function z'f^^ {y; k) for families of n vicious walkers, the i-th starting 
at (0, 2i — 2) and ending at {t,y + 2i — 2), i = 1,2, ... ,n, none of them running below 
the a;-axis. Our starting point is the double sum formula for zl^\y; n) in Theorem 8. 
It enables us to compute (in principle) full asymptotic expansions for all k > (see 
the Remark after the proof). The technique which we employ is singularity analy- 
sis. Our analysis reveals a phase transition at k = 2. We present the first term in 
the corresponding asymptotic expansions in Theorem 12 below. Previous asymptotic 
results for y = were found by Owczarek, Essam and Brak in [41, Eq. (4.58)], who 
gave predictions for the order of magnitude of the partition function Zg"^ (0; k) as r 
tends to oo. Our results confirm their predictions, but make them at the same time 
more precise, as we are also able to provide the multiplicative constants and, as we 
already mentioned, in fact, full asymptotic expansions. 

Theorem 12. Let y he a fixed non-negative integer. As t tends to oo, the partition 
function {y; n) is asymptotically 

2nt^-n(2n+l)/222n^-t+l^^-n/2 _ 1)! TT ^' (^ + ^ + ^^-1)! 

^ ^ -LI {y + 2iy. 



i=0 



X 



n fn\(y+h-2\ n 



it is 



2nt^-n(2n-l)/222n^-|n+l^-n/2 ^' + ^ + (l + 0{t-')) , if K = 2, 



.=0 (y+^^y- 

and it is 

nn-l^ \ * 9(n-l)(4n-5)/2^ _ o\2n-l 

^ ^ \ +-(n-l)(2n-l)/2^ H[H - 



(7.2) 



^/^^J 7r("-l)/2(^«-l)f+2n-l 



n-2 

X 



.=0 (^ + 2^)! 



Remark. This result shows different asymptotic behaviour for k < 2, for k = 
2, and for > 2. For k, < 2 the order of magnitude of the partition function is 
2nt^-n(2n+i)/2 growth rate is 2"', and the critical exponent is n{2n+ l)/2). 



31 



for K = 2 it is 2"*t-"(2"-i)/2, and for k > 2 it is (^^^^ t-("-i)(2n-i)/2 ^^le 

growth rate is not constant anymore, but grows with k), everything else is constants. 

Proof. We want to determine the asymptotic behaviour of (6.3). Clearly, for the 
product in front of the summation we just have to apply Stirling's formula. The result 
is that 

(^-1)! (t + 2z)!z!(y + n + i-l)! 



(l±jy +2n-2)! Al (t±^+n + ,_i)!(l^ + ,)!(y + 2z)! 

^ 2{n-l)t^-{{n+l){2n-l))/222n^-^n+^^{l-n)/2 

as t tends to oo. 

From now on we concentrate on the double sum in (6.3). Let us write t = 2r + y in 
the sequel. Since y is fixed and the summation index k comes from a bounded range, 
we may use the estimation 



ir + y)k 



C-i^ + n-k)^ {r + n-k)k 
to approximate the double sum by 

r n— 1 



l + 0(r-^) (7.5) 



5:i:MK"r)e^r^)f7-7-'r) 

z — k — 

■ i'L'-M"-') - * - 1" (1 + 

As is well-known, the most convenient method for determining the asymptotic 
behaviour of a sequence of numbers is singularity analysis as developed by Flajolet 
and Odlyzko [18] (see [19] for an introduction to that method). This method requires 
the generating function for the sequence of numbers to be known as a starting point. 
In our case, this means that we need a compact expression for 



r n—1 

k 



r=0 e=o k=0 



n - 1\ /y + 2k - 2\ l'2r + y - i - 1 
k )\ 2k )\ r-£-k 



We shall make use of the fact that 



r=0 
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where C{x) is the generating function (2.8) for Catalan numbers. This identity foUows 
readily from (2.7) by writing the left-hand side of (7.7) as 



°° ^2r + m^. ^ 
r 



r=0 

oo 

/ / '/ r -4- m \ I y r -A- m \ I y r -A- m \ I y r -I- m \ \ 

T 

r I \ r — 1 J \r — 1/ \ r — 2 

r=Q 

\m+l 1 ™/^/'™^m+3 1 \m+5 



y~v / /2r + m\ /2r + m\ /2r + m\ /2r + m\ ^ 



a;' 



= C(x)'"+^ + a;C(x)"^+^ + a;^C(a;)"^+^ + • • • . 

If we interchange summations in (7.6) so that the sum over r becomes the inner-most 
sum, and then apply (7.7), then we obtain the expression 

fn-l\fy + 2k-2\fy + i + 2n-2 



k J \ 2k J \ 2n-2k-l 



{2n -2k- 1)! m'—^^^^J^Iy^-'^'- 



Now we write 



^ + £ + 2n-2\ ^'"^ ^ U + h\(y + 2n-h-?,\ 
2n-2k-l J ^ \ h J\2n-2k-h-lJ ^ ' ^ 

(which follows from the Chu-Vandermonde summation formula (6.7)) and interchange 
summations so that the sum over £ becomes the inner-most sum. Using the binomial 
theorem and (2.10), we end up with the expression 

n-l 2n-2fc-l _ _ 2k - l\ fy + 2n - h - 3 



k=0 h=0 



k / \ h J \ 2n — h — 1 



Oh+l^ ^y+2kry( ^2\y+2k-l 

■hli2n-h-iy. ^ ^l^^ ^ ^ (7.9) 

'^^ VI - 4^2 (i + ^Vl - 4^2^ 

for the generating function (7.6). 

Singularity analysis now allows to extract the asymptotic behaviour of the coeffi- 
cient of 2;^'""'"^ as r tends to 00 out of the behaviour of the function in z at its singular 
points. We have to distinguish between three cases: k < 2, k = 2, and k > 2. 

Let first K < 2. Then, in the disk 1^1 < 1/2, the only singularity of the expression 
(7.9) is at 2; = 1/2. For z close to 1/2 we have the singular expansion 

g(^2)y+2fe -l ^^^^ 4^2)-l/22.+2/c-l +jj^o({l- 4z')'/^) , 

Vn^ (1 + 2^Vn^) 

(7.10) 
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where D is an explicit constant independent of z, whose value is however without rele- 
vance here. The coefficient of ,x"^ in (l-4x)-i/2 behaves like (4"^/v^)(l + 0(m-i)) 
as m — > oo. The transfer theorems [19, Theorems 5.4 and 5.5] then imply that the co- 
efficient of z^'^~'^^ in the left-hand side of (7.10) (which is what we need when we want 
to consider the coefficient of z'^^'+y in (7.9)) behaves like (22^+2^- V0rr)(l + 0{r-^)) 
as r — > oo. If we substitute this in (7.9), then we obtain that the coefficient of 2;^^+^ 
in (7.9) is 



n-l 2n-2k~l 

E E (-1)' 

k=0 h=0 



n — 1 
k 



2n - 2k - 1\ + 2n - h - 3 
h )[ 2n-h-l 



h\{2n-h-l)\ 



TTT (2 - K)f'+ 



as r tends to infinity. We interchange sums and write the (now) inner sum over k in 
hypergeometric notation. This gives 



2n-l 

E 

h=0 



2n-l\ fy + 2n - h-3 



h 



h\ {2n-h-l)\ 



2n-h-l 

22r+y+h^ 



7tt(2 - k)' 



—2F1 



2 + 2 n^l + ^ 



n 



n 



(l + 0(r-^)). 



Clearly, the 2-^1-series can be evaluated by means of the Chu-Vandermonde summa- 
tion (6.7). As it turns out, it is only non-zero for h > n — 1. After substitution of the 
result of the evaluation in the above expression, after replacement of /i by 2n — /i — 1, 
after taking into account that t = 2r + y, and finally combining the result with (7.4), 
the claimed result (7.1) follows upon little rearrangement. 

If K = 2, then it suffices to apply Stirling's formula to the closed form expression 
(6.10). 

Finally let k > 2. Then, in the disk 1^1 < 1/2, the "dominant" singularity (i.e., the 
singularity with least modulus; cf. [19]) is 2 = ^/k — 1/k. For z \/ k — 1/k we have 

y+2k-l 



C{z 



2\y+2k-l 



VT^^ (1 + ^vT^TiF) 



h+l 



K \ {K-2Y ) \ «-l ) 

(_l)/i+l«2/+2fc(^_ 2)2/1+1 
2/1+1 (k- l)2/+2fc+/. h _^2Ji!_j (7.11) 



The transfer theorems [19, Theorems 5.2, 5.4 and 5.5] imply that the coefficient of 



^2r-2fe 



in the left-hand side of (7.11) behaves like 



^y+2k(^ _ r,\2h+l 
( l)ft+l ^ ' 



K 



2h+l _ l^y+2k+h _ 1 
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Hence, the asymptotically dominating terms result from the summands in (7.9) where 
h is maximal, i.e., where h — 2n — 1, which, because of the binomial coefficient 
^2n-2fe-i^ occurring in (7.9), in turn forces k to be zero. If this is combined with 
(7.4), and if we again take into account that t — 2r + y, then the claimed result (7.3) 
follows upon little simplification. 

This completes the proof of the theorem. □ 

Remark. Since Stirling's formula (sec [52, Sec. 12.33]) does in fact provide a full 
asymptotic expansion for factorials, as well as do the transfer theorems [19, Theo- 
rems 5.2, 5.4 and 5.5] provide full asymptotic expansions, the above approach docs in 
fact allow to compute full asymptotic expansions for the partition function Z2^\y; k) 
in all three different regions, if needed. The reader should note that, if we want to 
compute more terms in the asymptotic expansion, we will have to extend the as- 
ymptotic expansion (7.5). The most convenient way to do that (for the subsequent 
generating function computations) is in the form 



C-i^+n-k)^ {r + n-k)k 

= ,^ Hy-n + k) ^ +0(r-"-^), (7.12) 

r + y + k {r + y + k)^ {r + y + k)^ ^ ^ ' 

with the Cj's appropriate constants independent of r. 

8. Exact and asymptotic results for the number of contacts of n vicious 
walkers. By definition, the mean number of contacts with the wall for families of n 
vicious walkers, the z-th starting at (0, Oj) and ending at (t, e^), none of them running 
below the x-axis (the wall), is 



Mi") (a - e; «) := 5^ N^^^ (a ^ e; £) £ k\ 
i>i 



where A^^*^*^^ (a — > e; £) is the number of these families of vicious walkers with exactly £ 
contacts with the wall. Clearly, we have 

Mi") (a ^ e; «) = «^^i"^(a - e; «), (8.1) 

with zj:""^ (a — > e; k) the partition function for these vicious walkers as defined in the 
Introduction. In addition, the normalized mean number of contacts with the wall is 
defined by 

— TT^ —■ 8-2 

Z^)(a^e;^) 

In this section we shall again concentrate on watermelon configurations which 
start on the wall and have deviation y, i.e., n vicious walkers as above with a = 
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(0, 2, 4, . . . , 2n — 2) and e = {y,y + 2,y + 4, . . . ,y + 2n — 2), and analyse their (nor- 
malized) mean number of contacts. Specializing (8.2) to the above choices of a and 
e, the quantity that we want to study is 

^ 8.3 

where, as before, Z^^^ (y; k) denotes the partition function for these watermelon con- 
figurations. Let us denote the normalized mean in (8.3) by M^''\y; k). Wc could use 
any of our formulas for the partition function zj:^\y; k) (such as (6.3) or (6.6)) to 
find explicit representations of M^^\y; k) as a quotient of two sums (double sums in 
the case of (6.3)). 

As was the case for the partition function, also these formulae for the mean number 
of contacts do not simplify in general, while they do for k = 1, and if y = also for 

K = 2. 

Theorem 13. For k.= 1, the normalized mean number of contacts, M^^\y\ 1), is 
equal to 

n{y + 2n + l){t - y) 



1 + 



{y + n){t + y + 4n) 



Proof. The value of the partition function at k = 1, Zf"'\y;l), was already 
determined in (6.9). For computing the numerator in (8.3), we make use of the formula 
(6.6). If this formula is differentiated with respect to k, and subsequently k is set 
equal to 1, then it is only two terms in the sum which survive. Some simplification 
then leads to the claimed formula. □ 

Whereas there does not seem a "nice" formula for the mean number of contacts 
for K = 2 for arbitrary y, there is one if y = 0. 

Theorem 14. For y — and k = 2, the normaUzed mean number of contacts, 
M^^\0;2), is equal to 



/2r+2n-2\ 
V r+n-1 ) 



-2(n- 1). 



Proof. Using the representation of the partition function Z2^\0; k) = Z^^ {k) as 
given in (4.1), the normalized mean number of contacts is 



r-1 



M^^ (0; 2) = 2 + . 

A comparison of (4.1) and (6.10) shows that we already evaluated (implicitly) the 
denominator of this expression. Hence, the evaluation of the numerator of (8.3) will 
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be accomplished once we have evaluated the sum in the numerator. We begin by 
converting it to hypergeometric notation, and obtain 



16n (r — l)r-i 



1 + 2n, 2 - r 
3-2r 



(r-1)! 

Next we reverse the order of summation and obtain 

2^+1 (2n)^_i 



;2 



2 — r, r 1 
2 - 2n - r ' 2 



(r - 2)! 

for the numerator of (8.3). We now apply the contiguous relation 



a, c 



; z 



c-1 
a-1 



a — l,c 
c-1 



; z 



a — c ^ 

+ T2Fl 

a — 1 



a — 1, c 



; z 



Thus, the 2^1-series becomes the sum of the two terms 



)r+l 



(2n), 



(r-1)! 



^F^ 



1 — r, r 1 
1 - 2n - r ' 2 



2'^+2n(2n) 



r-1 



(r-1)! 



Fi 



1 — r, r 1 
2 - 2n - r ' 2 



Both of the 2-^1-series can be summed by means of Bailey's summation formula (cf. 
[50, (1.7.1.8); Appendix (III.7)]) 



a,l — a 1 
b '2 



r(|)r(i + f) 



r(t + |)r(i-f + |) 

Some simplification then yields the claimed formula. □ 



(8.4) 



We now embark on determining the asymptotic behaviour of the normalized mean 
number of contacts. Predictions on the order of magnitude in the case of two walkers 
have been already made earlier by Brak, Essam and Owczarek in [8, Sec. 4.4.2]. 
The next theorem solves the problem for an arbitrary number of walkers, thereby 
confirming these predictions and, at the same time, making them more precise. Again, 
our method allows in fact to obtain full asymptotic expansions (see the Remark after 
the proof of Theorem 12). 

Theorem 15. As t tends to 00, the normalized mean number of contacts 
Mf^\y;K) is asymptotically 



l + C + 0{t-^), ifK<2, 



.5) 



where C is the quotient 



K 



1 



h=0 



en 



(2-k) 



2n-h+l 




h=0 



(?) (2 - 



2n-h 



1.6) 
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(8.7) 



it is 

2^-2"n (^^^ V^t-2n-y + 2 + 0{t-^/^) , if « = 2, 
and it is 

(^2+ 2(«-l)(«-2) + )' ^^'^>2- 

Remark. Thus, for the length of the walks being large, the normalized mean 
number of contacts is proportional to a constant if k < 2, it is proportional to the 
square root of the length of the walks if k = 2, and it is proportional to the length of 
the walks if k > 2. 

Proof. We start by observing that Theorem 8 implies that the normalized mean 
number of contacts Mj-^\y; k.) can be rewritten as 

^i"^(^;-) = i + S!r^' (8-8) 

where v}^\y; k) is the double sum in (6.3), and where 

{t-y)/2n-l 



e=o k=o ^ ^ ^ / \ 2 



k 



y + e+2n-2\ C-^), (2n - 2fe - 1)! ^1 ^ ^J+^ ^ 



2n-2k-l J (i_l+n-fc)^ 

We already determined the asymptotic behaviour of the denominator of the fraction 
in (8.8) in the proof of Theorem 12. What we still need is the asymptotic behaviour 
of the numerator. We attack this problem in the same way as the corresponding 
problem for the denominator. Here we will also have to use singularity analysis in the 
case that k = 2, because there is no closed form result available for the normalized 
mean number of contacts in that case (as opposed to for the partition function, which 
allowed us to do the asymptotics for k = 2 in Theorem 12 by just making use of 
Stirling's formula). 

Since for the 2 and > 2 we aim at computing an additional term (be- 

yond the leading term) in the asymptotic expansion, we will have to use the estimation 
(7.12), with N — 1, as a first step. After the substitution of this in the numerator 
of the fraction in (8.8), we compute the generating function with the coefficient of 
being the sum in the numerator. To be precise, again writing t = 2r + y, we want to 
find a compact expression for 



r n— 1 



r=0 £=0 k=0 



n - 1\ /y + 2k - 2\ /'2r + y - e - 1 
k )\ 2k [ r-i-k 



I A:=0 

l + i + 2 

2n-2k-l J^~'~ -ry—ry' ^ y _^ 



^ + ^ + 2" - 2\2n - 2/c - 1)! {2k)\ ( 1 + t^^^) (g.Q) 
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Now we perform the same steps as in the proof of Theorem 12: we interchange 
summations so that the sum over r becomes the inner-most sum, then sum over r 
using (7.7) (respectively (2.7) for the additional term), next use (7.8), then interchange 
summations so that the sum over i becomes the inner sum, and finally evaluate the 
remaining sum over £ by means of the binomial theorem. The result is that our 
generating function (8.9) turns into 



n-l 2n-2fc-l 

E E (-1)' 

k=Q h=0 



n - 1\ /2n -2k-l\ fy + 2n - h - 3 
k )\ h [ 2n-h-l 



{h + l)\{2n-h-l)\j,^-:^ 



plus a similar term resulting from the fraction k{y — n + k) / (r + y + k) occurring in 
(8.9). 

Now we apply singularity analysis to the function on the right-hand side. The 
considerations are entirely analogous to the corresponding ones in the proof of Theo- 
rem 12, except that here we have to consider one more term in the singular expansion. 
As it turns out in the end, the additional term k{y — n + k)/{r + y + k) does in fact 
produce smaller terms as we need. We leave the details to the reader. 

If, finally, the obtained results are combined with the corresponding ones for the 
sum in the denominator of the fraction on right-hand side of (8.8) which were found 
in the proof of Theorem 12, then after little simplification the claimed asymptotic 
expressions are obtained. □ 

Remark. The sums which appear in (8.5) in the case that k, < 2 can be expressed 
in terms of hypergeometric 2-P'i-series. (The sum in the denominator is a 2-^1-scries, 
while the sum in the numerator can be written as a sum of two 2-^1-series.) From 
there we see that they cannot be simplified, with the exception of a few special 
cases where 2-^1 -summation formulas are available (such as, for example, (6.7) and 
(8.4)). For instance, if k = 1, then the Chu-Vandermonde summation formula (6.7) 
is applicable. However, in that case it would have been easier to do the asymptotics 
directly from the closed form expression in Theorem 13. 
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